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Abstract 

In a previous paper [V. Delgado and J. G. Muga, Phys. Rev. A 56, 3425 
(1997)] we introduced a self-adjoint operator T(A) whose eigenstates can be used to 
define consistently a probability distribution of the time of arrival at a given spatial 
point. In the present work we show that the probability distribution previously 
proposed can be well understood on classical grounds in the sense that it is given 
by the expectation value of a certain positive definite operator J*^+^ (X) which is 
nothing but a straightforward quantum version of the modulus of the classical 
current. For quantum states highly localized in momentum space about a certain 
momentum pq ^ 0, the expectation value of J^"*"-* (X) becomes indistinguishable 
from the quantum probability current. This fact may provide a justification for 
the common practice of using the latter quantity as a probability distribution of 
arrival times. 



I. INTRODUCTION 

Standard quantum mechanics is mainly concerned with probability distributions of 
measurable quantities at a given instant of time. Such distributions can be inferred 
from the formalism in terms of projections of the instantaneous state vector \il^{t)) onto 
appropriate subspaces of the whole Hilbert space of physical states. However, one may 
also be interested in the probability that a certain physical quantity takes a definite 
value between the instants of time t and t + dt. Let us assume this quantity to be 
the position of a particle and restrict ourselves to one spatial dimension. What is the 
probability distribution of arrival times at a detector situated at a given point x = XI 
Standard quantum theory is based on the assumption that, by reducing sufficiently 



the experimental uncertainties involved, the outcomes of any measurement process 
will reproduce, within any desirable precision, the ideal distribution inferred from the 
spectral decomposition of a certain self-adjoint operator associated with the physical 
quantity under consideration. It is therefore implicitly assumed that the quantum 
formalism can provide a prediction for the experimental results without having to 
make reference to the specific properties of the measuring device involved. Since the 
distribution of arrival times at a given spatial point is, in principle, a measurable 
quantity that can be determined via a time-of-flight experiment, it is reasonable to ask 
for an apparatus-independent theoretical prediction. 

In classical statistical mechanics the above question has a definite answer: The 
(unnormalized) probability distribution of arrival times at X for a certain statistical 
ensemble of particles of mass m, moving along a well-defined spatial direction (i.e., 
either with momenta p > or with momenta p < 0) , is given by the average current at 

{J{X)) = fff{x,p,t)^ S{x - X) dx dp, (1) 
J J m 

where f{x,p,t) represents the phase-space distribution function characterizing the sta- 
tistical ensemble. In quantum mechanics, however, things turn out to be much more 
involved. In particular, a straightforward application of the correspondence principle 
would lead us to consider the expectation value of the current operator 

J{X) = ^(P\X){X\ + \X){X\P) (2) 

(P denoting the momentum operator) as the most natural quantum candidate for the 
probability distribution of the time of arrival at a point X. However, even though such 
a definition has been widely used in recent times [Q-Q, it cannot be considered as a 
satisfactory solution because of the fact that the expectation value of J{X) is not pos- 
itive definite, even for wave packets containing only positive-momentum components. 
Nonetheless, when quantum backflow contributions become negligible, one expects the 
expectation value of the current operator to be a good approximation to the actual 
probability distribution of arrival times. 

The difficulty for defining such probability distributions is nothing but a mere aspect 
of the more fundamental problem of the nonexistence of a quantum time operator 
conjugate to the Hamiltonian. The reason for this latter fact lies, basically, in the 
incompatibility of such a time operator with the semibounded nature of the Hamiltonian 
spectrum 

In spite of detailed work by Allcock Q denying the possibility of incorporating the 
time-of-arrival concept in the quantum framework, more recently there has been con- 
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siderable effort in defining a probability distribution of tfie time of arrival of a quantum 
particle at a given spatial point ||l|^ 0^121. The incorporation of such probability 



distributions in the formalism of quantum mechanics has both conceptual and practi- 
cal interest. In particular, this issue is closely related to the problem of the temporal 
characterization of tunneling (the so-called tunneling time problem) [p!3|-pl[, whose 
understanding is important for its possible application in semiconductor technology. 

In Ref. H] , Dumont and Marchioro, primarily concerned with tunneling-time distri- 
butions, proposed the probability current as a quantum definition of the (unnormalized) 
probability distribution of arrival times at a point sufficiently far to the right of a one- 
dimensional potential barrier. Leavens has shown that this result can also be derived 
within Bohm's trajectory interpretation of quantum mechanics by making the assump- 
tion that particles are not reflected back through the point X (i.e., {'4){t)\J {X)\il){t)) < 
does not occur for any t). On the other hand, Muga et al. 0] have provided an opera- 
tional justification of such a definition by simulating the detection of incoming particles 
by a destructive procedure. More recently, Grot et al. |^ have faced the problem from 
a somewhat different perspective. These authors construct a suitable self-adjoint op- 
erator in order to infer a probability distribution of arrival times from its spectral 
decomposition. More specifically, starting from the classical equations of motion for a 
particle moving freely in one spatial dimension and solving for the time, they arrive at 
the operator 



m 



as a natural candidate for the time of arrival of a quantum free particle at the spatial 
point X. In the above equation Xii{t) and Pnii) denote the position and momentum 
operators in the Heisenberg picture, respectively, and are related to the corresponding 
Schrodinger operators by 

Ojj(t) = e^-f^o*/^ 6 e-^^o*/^, (4) 

where O stands for X or P, and Hq = P'^/2m is the Hamiltonian of the free par- 
ticle. The operator (^) has the interest that it represents a quantum version of the 
corresponding classical expression t{X) = m[X — x(0)]/p(0) obtained by straightfor- 
ward application of the correspondence principle and the canonical quantization method 
||2^ , which states that classical equations remain formally valid in the quantum frame- 
work provided that one makes the substitution of Poisson brackets by commutators 
{ , } — > l/ih[ , ] and interprets the classical dynamical variables as self-adjoint op- 
erators in the Heisenberg picture. Of course, whenever the classical expression under 
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consideration contains products of dynamical variables having a nonvanishing Pois- 
son bracket (as is the case for position and momentum), the mere application of the 
canonical quantization method does not guarantee the unambiguous construction of 
the corresponding quantum quantity. In fact, a specific symmetrization or quantiza- 
tion rule has been explicitly chosen in Eq. (^). For instance, another possible quantum 
operator obtained from the same classical expression via a different symmetrization 



rule (which has been previously introduced by Aharonov and Bohm |23|) is 




-1 777 777 r ^ -1 \ 

'J Ph(0) Ph(0) L ' '\) 

Unfortunately, despite the fact that the above two operators are exactly what one 
would expect by virtue of the correspondence principle, none of them is self-adjoint. 
To circumvent this difficulty. Grot et al. proposed a modified time operator such that, 
when acting on states with no zero-momentum components, it leads to the same results 
as the operator (^) defined above. 

In a previous paper Q we followed a different route: Guided by the fact that, in 
general, a self-adjoint time operator conjugate to the Hamiltonian does not exist, we 
instead looked for a self-adjoint operator T[X) with dimensions of time, conjugate to 
a conveniently defined energy operator having a nonbounded spectrum. We showed 
that the orthogonal spectral decomposition of such an operator can be used to define 
consistently a probability distribution of arrival times at a given spatial point within 
the standard formalism of quantum mechanics. 

In this paper we are mainly interested in relating the formulation proposed in Ref. 
with the corresponding classical formulation. A quantum expression looks more 
natural as long as it is possible to derive it from a known classical quantity by applying 
certain specific quantization rules (even though such a procedure is by no means a 
necessary condition for the validity of a quantum formulation). We begin by briefly 
reviewing the relevant formalism in Sec. II. In Sec. Ill we consider the semiclassical 
limit of the proposed probability distribution of arrival times. In Sec. IV we show that 
such a probability distribution can be well understood on classical grounds in the sense 
that it is formally analogous to its corresponding classical counterpart. In this section 
we also provide a relation between the expectation value of the self-adjoint operator 
T(X) and the expectation values of the operators T[X) given by Eqs. (^) and (^). 
In Sec. V we analyze under what circumstances the proposed probability distribution 
can be replaced, to a good approximation, by the probability current, which has been 
frequently used, in practice, as a quantum probability distribution of arrival times. 
Finally, we conclude in Sec. VI. 
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II. FORMALISM 



In looking for a probability distribution of arrival times within the framework of 
standard quantum mechanics, we introduced in Ref. a self-adjoint energy operator 
TL defined by 

n = sgn(P) (6) 
Im 

where, in terms of a basis of momentum eigenstates {p)}, the operator sgn(P) reads 

sgn(P)=/ dp{\p){p\-\-p){-p\). (7) 





The normalization has been chosen so that the states \p) satisfy the closure and or- 
thonormalization relations 

+ 00 

dp\p){p\ = l, (8) 

-oo 

{p\p') = 5{p-p'). (9) 

The motivation for introducing the operator 7i, which essentially represents the 
energy of the free particle with the sign of its momentum, lies in the fact that, unlike 
the Hamiltonian, it exhibits a nonbounded spectrum. It is therefore possible to define 
a self-adjoint operator with dimensions of time 'T{X) by simply demanding it to be 
conjugate to TL, i.e., 

[n, f{x)\ = e-'^^/^ [n, t(o)] e+'^^/^' = in. (lo) 

This procedure led us to an operator T{X) whose orthogonal spectral decomposition 
reads 



T{X)= dTr\T;X){r;X\, (11) 

J —oo 

where the states \t;X), which constitute a complete and orthogonal set, are given by 



t; 



m 



In order to facilitate an interpretation in terms of measurement results, it turns 
out to be most convenient to decompose the eigenstates |t;X) as a superposition of 
negative- and positive- momentum contributions, in the form 
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\t-X) = \t = -T,--X) + \t = +T,+]X), (13) 
with |t, ib;X) defined by 

%±-X) = h-'/^ rdp,/^e^(i^*^''^)/^ \±p). (14) 
Jo \ m 

As can be easily verified, even though the states |t, ±; X) constitute a complete set they 
are not orthogonal. Specifically, 

/ + CO 
dt \t,a-X){t,a;X\ = 1, (15) 

(t, a|t', a') = |(5(t - t') - P^^^^^l . (16) 



Despite this fact, the decomposition (13) is interesting because the variable i, unlike 
r, admits a proper interpretation as a physical time. In particular, the states |t,ib;X) 
not only have the desirable time-translation property 



JHot'/h 



t,±-X) = \t + t',±;X), (17) 



but also transform under time reversal as |t,ib)— >|— t,=p). 

Consider a free particle propagating along the x axis toward a detector located at a 
given point X. We shall assume that its actual state at t = is, in the position repre- 
sentation, either a linear superposition of positive plane waves |^/^+(0)) (corresponding 
to particles arriving at the detector from the left) or a linear superposition of negative 
plane waves \ip-{0)) (corresponding to particles arriving at the detector from the right). 
At any instant of time the state vectors |V'±(^)) satisfy the identity 

|V^±(t))^G(±P)|V'±(t)), (18) 

where 0(ibP) are projectors onto the subspaces spanned by plane waves with posi- 
tive/negative momenta 



G(±P) = / dp \±p){±p\ . (19) 

and vanishing (in 



It can be shown that for normalizable states satisfying Eq. ([ 
momentum representation) faster than p as p — > 0, it holds that [|/ 
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dTr(V± |±t;X)(±t;X| V±) 

-oo 

I^::drT{i;^iT)\J{X)\ij^iT)) 



-oo 

/^dr(V^±(r)|J(X)|^±(r)) 



(20) 



where use has been made of Eq. (|ri|) and \ip±) denotes the state of the particle in the 
Heisenberg picture, 

l^±> = e*'^°^/"|V±(T)) = |V±(0)). (21) 

As already stated, the right-hand side of Eq. (^) can be recognized as a quantum 
version of the mean arrival time at X, obtained by straightforward application of the 
correspondence principle to the analogous classical expression for a statistical ensem- 
ble of particles propagating along a well-defined spatial direction. Furthermore, the 
positive-definite quantity {^jJ± \±t; X){±t; X \ ip±) satisfies 

/ + 00 
dT{^±\±T■,X){±T■,X\^lJ±) = (V'ilV'i) = 1- (22) 
-oo 

Therefore, in a quantum framework, the mean arrival time at X can be defined consis- 
tently by 

/ + 00 
dTT{^l^±\±T■,X){±T■,X\^l;±). (23) 
-oo 

Accordingly, the probability amplitude ^±{t = t;X) of arriving at X at the instant 
t = T, coming from the left/right, would be given by 



^±{t = T;X) = {±T■,X\^P±) = {t = ±T,±;X\i;±) 

= /.-V2rdp#(±p|^±)e-<^^H/^ (24) 
JO V m, 

and the corresponding probability density takes the form 
Jo JO mh 

(25) 

The above formulation can be generalized in order to define a probability distri- 
bution of arrival times at an asymptotic point X behind a one-dimensional potential 
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barrier. Indeed, provided that the potential V{x) vanishes sufficiently fast, far away 
from the scattering center, as to guarantee the validity of the standard scattering for- 
malism, it can be shown that Eqs. (p2D-(p5[) remain formally valid with the only 
substitution 

where the (unnormalized) freely evolving transmitted state IV'tr) can be written in terms 
of the scattering operator S as 

POO 

|^tr> = e(P)S|Vin) = / dpT{p){p\i^i^)\p). (27) 
Jo 

In the above equation T{p) denotes the transmission coefficient characterizing the po- 
tential barrier, and the state vector iV'in) [which is assumed to satisfy the identity 
iV'in) = 0(-P) iV'in)] represents the incoming asymptote of the actual scattering state 
of the particle at t = 0. Whenever this latter state |^(0)) does not overlap apprecia- 
bly with the potential barrier, it becomes physically indistinguishable from and, 
consequently, it is not necessary to discriminate between them in practice [p^ . 

Since the presence of a potential barrier is not relevant for our purposes in this work, 
in what follows we shall restrict ourselves to the free case. More specifically, we shall 
consider a freely moving particle characterized by a state vector |V'±(i)) satisfying Eq. 
(^). Nonetheless, this assumption does not imply any loss of generality in practice since 
the formulation below can be systematically generalized by means of the substitution 



III. SEMICLASSICAL LIMIT 



In spite of the fact that the expectation value 

± (J(X))± ^ ±(^±(r)|J(X)|V±(r)) = ±{^l;±\Jn{X,T)\ij±) (28) 

cannot be properly considered as a probability density of arrival times, it represents, 
however, a natural quantum version of the corresponding classical probability density. 
For this reason it is instructive to investigate the connection between such an expec- 
tation value and the quantity {ip± \±t;X) {±t; X\7p±) , which, as stated above, can be 
interpreted consistently as a quantum probability density of arrival times. To this end, 
by inserting twice the resolution of unity in terms of a momentum basis, we write 
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/^eTi(p-p')^/?i. 



(29) 



Expressing next the probability amplitude {^p\tjj±) in polar form as 

(±p|V^±) = |(±p|V'±)|e*'^±(^')/'* (30) 
and introducing the functional 

poo 

I± [f] ^ / dp/(p) |(±p|V±)| e*^±(P)/^ (31) 
Jo 

where the phase x±ip) is defined as 

2 

X±{p)^(l>±ip)-^r±pX, (32) 
we can finally rewrite Eq. (p9|) in the form 



± (V'il Jh(X, t)|^±) = ^1{I*± \P] I± [1] + c.c.) . (33) 



On the other hand, the probability density (|25|) (which, unlike the probability current, 
is manifestly positive definite) can be written in a completely analogous manner as 

(V^± |±r; X) (±r; XjV'i) = ^ (4 [Vp] ^± [^^) • (34) 



The two expressions ( p3|) and ( |3^ ) are especially suitable for investigating the semi- 
classical limit ?i ^ 0. Indeed, in this limit the asymptotic expansion of the Fourier-type 
integral ( |3l| ) is given, to leading order, by |25| 



I± [/] ~ ^s'^[^±(^'°)l e*^±(«')/V(Po) KipolV'i)! Jr^J^: (35) 

V lx±(Po)| 

where po is the stationary point of the phase x±{p)i defined implicitly by 

x'±{po) ^ Mpo) -^t±X = 0, (36) 
ni 

and primes are used to denote differentiation with respect to the momentum p. 



By substituting Eq. (^) into Eqs. (33) and ( p^ ) we obtain, to leading order as 
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(V'i \±t;X) (±r;X|V^±) ~ ±(^±| Jh(X, r)|V^±) ~ t)|2. (37) 

m 

In the last step of the above formula we have used that 

V^±(X,t) = (X|eC-A|V;±) = h-^/^I^ [1] , (38) 
so that, as ?i — > we have 

|V'±(X,r)| . (39) 

lx±(Po)| 

The information contained in Eq. ( pTj ) represents the main result of this sec- 
tion. This equation reflects that the proposed probability density of arrival times 
Kibr; Xj-;/;-!-)!^ coincides, in the semiclassical limit, with the quantum probability cur- 
rent ib(^/;-|-| Jh(X, r)|'i/'±), which in turn is given in this limit by the product of the 
probability density \ip±{X,T)\'^ and the group velocity Po/m (and, consequently, be- 
comes a positive quantity). This fact suggests that the probability density defined in 
Eq. (|2^) is nothing but a quantum version of the modulus of the classical average 
current {J{X)), which, as already stated, plays the role of a probability distribution 
of arrival times at X for a classical statistical ensemble of particles moving along a 
well-defined spatial direction. In the next section we shall see that this is indeed the 
case. 



IV. POSITIVE-DEFINITE CURRENT 



Let us concentrate on Eq. (12), which defines the eigenstates of the operator T{X) 
in terms of a basis of momentum eigenstates. By introducing the self-adjoint operator 



/-1-00 , 
dpJ\p 
-oo 



(40) 



we can express |r; X) in an alternative form that exhibits no explicit dependence on any 
particular representation and proves to be most convenient for our purposes. Indeed, 



\r;X) 



dp • 



\P\ 



m 



\P\ 



m 



(41) 
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Correspondingly, the probability amplitude for particle detection at the spatial point 
X [coming from the left(+)/right(— )] at time t = T takes the form [Eq. (24)] 




^±{t = r-X) = {X\ \li:-le-'^-'^\i:^) = (X| ^ ^ |V^±(t)). (42) 
This equation shows that the probability amplitude of arriving at X at time r is nothing 



but the probability amplitude of finding the state \j\P\/m\X) in the (Schrodinger) 
state vector \'4^±{t)) characterizing the particle dynamics at t = r. On the other hand, 
the corresponding probability density reads 



P^x\r) , = \^±{t = T,Xr = (V'iWI \I^S{X -X) W^|^±(r)), (43) 

quant \ m \ m 

where use has been made of the identity = 6{X — X). 

Before proceeding further it is convenient to consider the normalized probability 
distribution of arrival times at the point X for a classical statistical ensemble of free 
particles of mass m, coming either from the left [p > 0) or from the right {p < 0). Such 
a probability distribution can be obtained from Eq. (||). Indeed, by defining J^~^\X) 
as the modulus of the classical current 

j(+)(X) = \ J{X)\ = \p\/m6{x - X) (44) 
and using a convenient notation, it takes the form 




\P\ 



, =±(J(X)>± = /'//±(i-,p,T)Ma(i_X)<!i<ip = {j(+>(X))±, (45) 

class J J m 

where the phase space distribution function satisfies the identity 

f±{x,p,t) = e{±p)f±{x,p,t) (46) 

and the modulus in the integrand of Eq. ( [45| ) comes from the normalization factor 
(which, in the free case, takes the value ±1). 

A comparison between expressions (^) and (|4^) shows that the quantum proba- 
bility density of arrival times defined above [Eq. (^)] can be considered as a quantum 
version of the corresponding classical expression, obtained by associating to the average 
of the classical positive current J^^\X) = \J{X)\ the expectation value of the positive 
definite current operator J^'^\X) 
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j(+) (^x) = ^-^ 6{x - X) ^ J(+) {X) = \ — 6{X -X)\r-^. (47) 




m \ m 

It should be noted, however, that the relation existing between the classical current 
J{X) and the corresponding quantum operator J{X) [given by Eq. (^)] is somehow 
different from that existing between J^~^\X) and J^~^\X). Indeed, J{X) can also be 
considered as the quantum operator corresponding to the classical current J{X) by 
virtue of the Weyl-Wigner quantization rule, whereas the same does not hold true for 
the positive current defined above. 

The Weyl-Wigner quantization rule is a mapping that associates with every phase- 
space function g{x,p) a quantum operator G{X, P) with an expectation value satisfying 



= JjMx,p)g{x,p)dxdp, (48) 
where the Wigner function f-\j\i{x,p) plays the role of a quasiprobability distribution 



function in phase space |26| and can be expressed in terms of the quantum density 



operator p characterizing the physical system as [27| 



By substituting Eq. ( |49| ) into Eq. (^) it can be shown that G[X, P) is given by 

G{X, P) = ^JJJJ g{x, p)e'^0(X-g)+r{P~p)] (gp) 

and taking g{x,p) = p/m5{x — X) in the integrand of Eq. (pO|) one arrives, after 
some algebra, at the current operator J{X) defined by Eq. (||). Even though a similar 
relation does not exist for J^^\X) = \J{X)\ it still holds true that the positive-definite 
operator J^~^^ {X) represents a natural quantum version of the modulus of the classical 
current. Accordingly, for free particles propagating along a well-defined spatial direction 
the probability density of the time of arrival at a given point X at time t = t can be 
defined consistently, within both a classical and a quantum-mechanical framework, as 
the instantaneous mean value of the modulus of the current 



class 



(J(+)(X))±, (51) 



Prir) ^ = (V± \±t;X) {±T;X\i;±) = (V'±(r)| j(+)(X)|V'±(r)). (52) 

quant 



We shall next concentrate on the operator T{X). From the definition (|llj) one finds, 
taking Eq. ( pl| ) into account, that T(X) satisfies 
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± Q{±p)T{x) e{±p) = e(±p) 



+ C)0 



dTTj^\X,T 



where J^^^ {X, r) denotes the positive current in the Heisenberg picture 

jW(X,t) = e^^«^/^ J^+\X) e-^^"^/'*. 
Accordingly, the mean arrival time at X [ Eq. (p3|)] can be expressed as 



(53) 



(54) 



{tx)± = ±{i^±\r{x)\i;±) 



dTr(V'±(T)|j(+)(X)|V^±(r)). 



(55) 



The above equation gives the mean arrival time in a form that can be recognized as a 
quantum version of its classical counterpart in terms of the probability distribution (45). 
Indeed, the positive current {il^±iT)\J^^\X)\'tlj±{T)) enters Eq. (|55| ) as a probability 
density of the time of arrival at X. It should be stressed that contrary to what happens 
with Eq. (pO|) , by virtue of Eq. (|53|) the above equation is valid for any iV'it (''")) 
satisfying the condition (|l8|) . When in addition to this latter condition it is also satisfied 
that 



lirn (piV'zb) 
p— >±oo 



0, lini p ^ {p\tlJ±) 



0, 



(56) 



it can be shown, after some algebra, that the mean arrival time {tx)± can also be 
expressed in the alternative forms 



{t 



x)± 



{tx)± 

, 1 



+ 00 



drT(^±(r)|J(X)|7/.±(r)), 



(57) 



X-Xn{0) 



m 



+ 



m 



{t 



m 



Pr{0) Ph(0) 
^-^h(0)1 J 



X-XuiO) \iJ±), (58) 



m 



(59) 



In the derivation of the above formulas use has been made of the fact that the position 
operator transforms under spatial translations as 



^-^Px/n ^+^px/^ ^ -X) 



(60) 



and that, at t = 0, quantum operators in the Schrodinger picture become indistinguish- 
able from those in the Heisenberg picture, so that, in particular, we have X = Xh(0) 
and P = Ph(0). 
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Equations (^8|) and ( |59| ) give the mean arrival time at X in terms of the operators 
T{X) introduced by Grot et al. and by Aharonov and Bohm [Eqs. (^) and (|5|), re- 
spectively]. A connection between these operators and the self-adjoint "time" operator 
T(X) can be derived from a comparison with Eq. (|55|). Indeed, as long as conditions 
(p^) and ( |56[ ) are satisfied, the expectation value of 'T{X) coincides with the expecta- 
tion values of the operators T(X), which have the interest that they can be obtained 
by quantizing the classical expression t(X) = m[X — x{0)]/p(0) according to different 
standard quantization (ordering) rules. 

Under the same conditions, the expression of {tx)± given by Eq. (^) becomes 
also indistinguishable from that given by Eq. (|5^), which involves the usual quantum 
probability current and has been frequently used, in practice, as a quantum definition 
for the mean arrival time. This fact may provide additional justification for the latter 
expression, whose validity in a quantum framework might, in principle, be questionable. 
Indeed, despite the formal analogy between Eqs. ( |55|) and (|57|), they have a somewhat 
different physical meaning: While the positive current {X) is a positive-definite op- 
erator and its expectation value enters Eq. ( [55| ) playing the role of a probability density, 
the expectation value of the usual probability current J{X) can take negative values 
and consequently cannot be interpreted as a probability distribution of arrival times. 
Since usually this has been the case, however, it is instructive analyzing under what 
circumstances the expectation values of J{X) and J^'^\X) become indistinguishable. 
This will be the aim of the next section. 



V. PROBABILITY CURRENT VERSUS POSITIVE CURRENT 



In the preceding section we have seen that for state vectors satisfying Eqs. (|ll 
and (|56|) the mean arrival time {tx)± can be equally calculated by using the positive 
current or the probability current [Eqs. ( |55|) and (^), respectively]. Furthermore, by 



usmg 



dr e 



6{p -p) + -r-rS{p' +p) 
p\ \p\ 



in Eqs. (Pq) and (29), we obtain 

X) 

dr (r) I J(+) (X) I V± (r)) = ± / dr {^^ (r) \J{X)\^^ (r)) , 



(61) 



(62) 
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so that the total probabihty of arriving at X (at any instant of time) can also be 
calculated in terms of J^~^^{X) or J{X). Despite the interchangeable role that these 
quantities play in the above integral expressions, when the state vector describing the 
quantum particle exhibits a considerable spread in momentum space the contribution of 
quantum interference effects may be important and the expectation values of J{X) and 
J^^\X) can be appreciably different. This is most easily seen by considering a super- 
position of two nonoverlapping wave packets with well-defined momentum. Specifically, 
we shall consider a state vector given by 

1^/^+) = ai iV^i) + a2 IV'2), (63) 

where the coefficients ai, ai are real and (pIV'j) (j = 1) 2) are assumed to be minimum 
Gaussian wave packets centered at the spatial point xq, with momentum spread Ap 
and average momentum pj , respectively, 



27r(Ap)^ 



-1/4 



exp 



2Ap 



■PXo 



(64) 



We take p2 > pi > and Ap <^ {p2 — pi) in order to guarantee that the above two 
wave packets do not overlap appreciably. Under these assumptions, the probability cur- 
rent ('i/'+l Jh(-^5 ''")|V'+) and the probability density {iI)+\J^\x,t)\iIjj^) can be written, 
respectively, as 



(^+1 Jh(X, r)|V+) = ^\{I*\P]I [1] + c.c.) 



(V+|4+^(X,t)|^^ 



mh 



(65) 
(66) 



where now the functional / [/] is given by / [/] = Ii [/] + I2 [/] with 



2Tr{Apy 



-1/4 



dpf{p) exp 



P-Pj 
2Ap 



.pxo 



exp 



p"^ T ,pX 
2m h h 



(67) 

To obtain an analytical estimation (as a function of r) for the probability current 
and for the probability density of arrival times [as given by Eqs. (|6^) and (^), re- 
spectively] we shall next consider the asymptotic expansion of the above integral. To 
leading order as Ap — > 0, we have [25 



~ aj V2 



27r(Ap)2l'^V(Pi)exp 



-|Lj-zp,(xo-X)/n 



O 



(ApY 



(68) 
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By substituting Eq. 
asymptotic expressions 



into Eqs. (pBD and ([6^), one obtains, after some algebra, the 



2V27r 
mh 



Ap (ajpi + alp2 + 0102(^1 + P2) 



X cos (pI - p\)T/2m + {p2 - Pi)(xo - X) /fi + O (Ap) 



(69) 



(V^+|JW(X,r)|V'^ 



2V2vr 
mh 



Ap a^pi + a2P2 + OLxa2yJp\P2 



X cos 



\pI - pl)T/2m + ip2 - pi){xQ - X)] /h + [{Apf\ ) . 



(70) 



From Eq. ( pOD we see that the probabihty current can take negative values whenever 
the interference term dominates over both the first and the second one. It is not hard 
to see that this is the case when it holds that 



02 pi 



(71) 



and, under these circumstances, Eq. (^) is given, to a good approximation, by 



{^+\MX,t)\^, 



2V27r 
mh 



Apaia2P2Cos p\t /2m + p2{xo - X) /h + {Ap) 



(72) 

In contrast, the probability density (^4.| J^^'*(X, r)|^/'_|_) remains always positive. This 
fact, which is evident from Eq. (|66[), can also be verified from the asymptotic expression 
(^) by noting that the dominance of the interference term would require that 



1 « —J'— <. 1, (73) 

02 V P2 

which is obviously impossible. 

On the other hand, it can be readily verified that for nonnormalizable states with 
a well-defined momentum = \p) {p ^ 0), the quantity {il^±{T)\J^'^\X)\ilj±{T)) 

becomes indistinguishable from the usual probability current ziz{ip±{T)\J{X)\ip±{T)) . 
One expects this fact to be also true for normalizable states describing particles with 
a highly defined momentum, and this is indeed the case as can be inferred from the 
asymptotic behavior (as the momentum uncertainty approaches zero) of the integrals 
defining the expectation values of J{X) and J^~^\X) [see Eqs. (|69|) and (|70|) above and 
take Q2 = 0]. 
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In scattering problems one is usually concerned with particles propagating with a 
well-defined velocity toward a localized interaction center. Such particles are character- 
ized by quantum states highly concentrated in momentum space about a certain mo- 
mentum po 7^ 0. Under these conditions, the expectation values of J{X) and J^'^^X) 
coincide to a good approximation, so that the probability current yields essentially 
correct results for the probability density of the time of arrival at a given point. This 
fact provides a justification for the use of J{X) in this kind of problems. 

Finally, it is worth noting that Eq. (57), which gives the mean arrival time {tx)± 
in terms of the probability current, is applicable even for state vectors having a large 
momentum uncertainty. Indeed, its validity only requires the fulfillment of conditions 
(H) and (H). 



VI. CONCLUSION 

Quantum theories are assumed to be more fundamental in nature than the corre- 
sponding classical theories. Consequently, it is possible, in principle, to define quantum 
quantities without resorting to the correspondence principle. In practice, however, 
the correspondence principle proves to be extremely useful in the construction of the 
quantum counterpart of a certain classical quantity. For instance, the canonical quan- 
tization method represents an invaluable tool for the construction of quantum field 
theories. Furthermore, one usually gets a better understanding of a quantum theory 
when a clear and unambiguous relationship can be established between quantum and 
classical quantities. 

In this paper we have been particularly interested in investigating the connection 
between the expressions previously proposed in Ref. for the probability distribu- 
tion of the time of arrival at a given spatial point and their corresponding classical 
counterparts. In particular, we have shown that, in the semiclassical limit ?i ^ 0, 
the proposed probability density of arrival times coincides with the modulus of the 
quantum probability current. Indeed, Eq. ( |37| ) can be rewritten in the form 

K±r;X|V^±)|2~K^±(r)|J(X)|V^±(r))|. (74) 

This result has the interest that, at a classical level, the current of a statistical ensem- 
ble of particles propagating along a well-defined spatial direction plays the role of a 
probability distribution of arrival times. Therefore, Eq. ( |74|) reflects that the quan- 
tity |(ibr;X|?/;-|-)p has the correct semiclassical limit and suggests that it represents a 
quantum version of the corresponding classical expression. We have explicitly shown 
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that this is the case by expressing the probabihty distribution |(ibr;X|V'±)P as the 
expectation value of a certain positive definite current operator. Indeed, by making use 
of the fact that the probabihty amphtude of arriving at X at time r coincides with the 



probabihty amphtude of finding the state vector y\P\/m \ X) in the Schrodinger state 
of the particle \iP±{t)), one can write the corresponding probability density in the form 



^ |(±T;X|fc)|= = (fc(T)|J<+>(X)|^fa(T)), (75) 

quant 



where the positive-definite operator 




jW(X) = y^<^(X-X)^/^, (76) 

can be immediately recognized as a straightforward quantum version of the modulus 
of the classical current |J(X)| = \p\/m5{x — X). The existence of a remarkable for- 
mal analogy between the corresponding classical and quantum expressions is therefore 
apparent: For particles propagating along a well-defined spatial direction, the proba- 
bility distribution of the time of arrival at a given point can be inferred, within both a 
classical and a quantum framework, from the mean value of the modulus of the current. 
On the other hand, for normalizable states satisfying the identity 

|V±(t))^e(±p)|V'±(t)) (77) 

and vanishing faster than p asp approaches zero, the mean arrival time at X, {tx)±, can 
be equally calculated in terms of the positive-definite current (X) or in terms of the 
standard probability current. This interchangeable role is not restricted to the mean 
arrival time. Indeed, we have seen that for physical states with a sufficiently well-defined 
momentum the expectation values of J{X) and J^'^\X) become indistinguishable, so 
that the probability current yields essentially correct results for the probability density 
of the time of arrival at a given point. This fact may provide a justification for the 
common practice of using the expectation value of J{X) in this kind of problem. 

Furthermore, under the same above conditions, the expectation value of the self- 
adjoint "time" operator T(X) coincides with the expectation values of the operators 
T{X) introduced by Grot et al. and by Aharonov and Bohm [Eqs. (|3|) and (^), 
respectively], which have the interest that they can be considered as straightforward 
quantum versions of the classical expression t{X) = m[X — x(0)]/p(0). 

In summary, we have shown that the formalism developed in Ref. |Q] for the time 
of arrival of a quantum particle at a given spatial point can be reformulated in a form 
that exhibits a remarkable formal analogy with the corresponding classical formulation. 
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Note added in proof: 

As stated in the Introduction, tlie idea that according to standard quantum me- 
chanics measuring results of physical quantities can be inferred from the spectral de- 
composition of a certain self-adjoint operator, without having to make reference to 
the specific properties of the measuring device involved, plays a central role in our 
treatment. In this regard it should be mentioned that a completely different view is 
developed in a recent preprint by Aharonov et al. |^^. These authors, by explic- 
itly modelling the measuring device, arrive at the conclusion that the time of arrival 
cannot be precisely defined and measured in quantum mechanics. In the view of the 
present author, however, this pessimistic conclusion can be mitigated, in part, by the 
assumptions on which it is based p9| . 

On the other hand, in connection with Eq. ( |7^ ) [which is only valid in the semi- 
classical limit] it is interesting to note that, as shown by McKinnon and Leavens j^, 
within the framework of Bohmian mechanics the modulus of the quantum probability 
current provides a definition for the probability density of arrival times which is of 
general applicability |^^. Therefore, while according to Eq. ([75|) within conventional 
quantum mechanics the probability density of the time of arrival can be defined as the 
mean value of the modulus of the current, within Bohmian mechanics it is the modulus 
of the mean value of the current the relevant quantity (the range of applicability of the 
two expressions is different, however). 
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